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Abstract 

The three-state Potts model with antiferromagnetic nearest-neighbor (n.n.) 

and ferromagnetic next-nearest-neighbor (n.n.n) interaction is investigated 

within a mean-field theory. We find that the phase-diagram contains two 

kind of ordered phases, so-called BSS phase and PSS phase, separated by a 

discontinuous phase transition line. Order-disorder transition is continuous 

for the weak n.n.n. interaction and becomes discontinuous transition when the 

n.n.n. interaction is increased. We show that the multicritical point where the 

order-disorder transition becomes discontinuous is indeed a tricritical point. 
PACS numbers: 64.40.Kw, 75.10.Hk 
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I. INTRODUCTION 



Antiferromagnetic three-state Potts model has interesting properties. It is described by 
the following Hamiltonian: 

H = J £ S(s it a s ) (1) 

where < ij > indicates summation over nearest neighbor pairs, and Sj = 1, 2, 3 denotes 
three-state Potts spin on the z'th site. In antiferromagnetic case (J > 0), the neighboring 
two spins prefer to take different values. Typical ground state configuration of the model on 
the square lattice is depicted in Fig.(l). One can change the state of certain spin in Fig.(l) 
from "2" to "3", or vice versa, without any energy cost. One half of all the spins are such 
"semi-free" spins, therefore ground state is infinitely degenerated. It should be noted that if 
we divide the lattice in Fig.(l) into two interpenetrating sublattices, there are only "1" spins 
on one sublattice, while random mixture of "2" spins and "3" spins is present on another 
sublattice. 

For the model on simple cubic lattice, it is known that long-range order is realized at 
finite temperature in spite of the high degeneracy of ground-state which usually leads 
to the disordered state at all temperature range ||. Recently this model (on the simple 
cubic lattice) has been studied intensively concerning two interests; one on the nature of the 
order-disorder transition and another on the nature of ordered phase. 

Order-disorder transition is continuous and the associated critical exponents are esti- 
mated by various methods, such as Monte Carlo simulation and the coherent- anomaly 
method @. 

As for the nature of ordered phase, it is known that so-called broken-sublattice-symmetry 
(BSS) phase is realized at a sufficiently low temperature in which one of the sublattices 



is dominated by one spin state, while the second sublattice is dominated by the mixture of 
the remaining two spin states. 

Recently several different claims have been made about the nature of the ordered phase at 
the region just below the transition temperature. Lapinskas and Rosengren concluded that 



so-called permutationally-symmetric-sublattice (PSS) phase, in which each two sublattices 
are dominated by one spin state, is realized at a very narrow temperature range just below 
the transition point, based on cluster- variation method analysis [|TT|,|i"2f and Monte Carlo 



simulation |T3" |. On the other hand, Kolesik and Suzuki have performed Monte Carlo simu- 
lations and observed a rotationally symmetric phase at a certain range near the transition 
point, finding no evidence of the PSS phase fl4[ . 



Now we consider the effect of next-nearest-neighbor (n.n.n.) ferromagnetic interaction. 
Let us consider the following Hamiltonian: 

H = J Yl S( Si , Sj )-<yJ £ 6{ Si , Sj ) (2) 

<i,3> 11 II <i,3> n n ?! 

where the first and the second summation runs over all nearest neighbor and next-nearest 
neighbor pairs on the simple cubic lattice, respectively. We assume that J > and 7 > 0. 
The lattice consists of two sublattices, which we refer as A and B. The first term in (0) 
represents the inter-sublattice antiferromagnetic interaction and the second one represents 
the intra-sublattice ferromagnetic interaction, which resolves the high degeneracy of ground 
state. 

Thus the n.n.n. interaction affects the nature of the ordered phase because the BSS phase 
costs energy proportional to 7J compared to the PSS phase. This effect for the models on 
the square lattice has been studied by several methods fl5HT7jl . 

In three dimensions, we expect another effect produced by a strong n.n.n. interaction; 
one can consider 7 — > 00 limit as two independent systems of the ferromagnetic q = 3 Potts 



model, which undergoes the first-order phase transition in three dimensions ||18|| . Thus the 
order-disorder transition is expected to become discontinuous as 7 becomes large. 

Mean-field theory gives qualitatively correct answer for both 7 = 0(AF) and 7 — > 00 (F) 
cases in three or greater dimensions, namely the BSS phase is realized below the continuous 
phase transition point in the AF case, and discontinuous phase transition occurs in the 
F case. We investigate the intermediate region using mean-field theory in the following 
sections. 



II. MEAN-FIELD ANALYSIS 



To perform a mean-field calculation, we use a method of variational free energy which 
is equivalent to equation of self-consistency [|19|. Let us consider the following mean- field 
Hamiltonian: 

3 

ff o = EEM s i) s ) ( 3 ) 

i s=l 

where hi s is a mean-field acting on the i'th spin. Boltzmann probability factor for the 
system described by H is denoted by P ({ s i})> which is a product of probability factors of 
individual spins: 

P ({ Sl }) = IkpM) (4) 

p .( s ) = exp(-(3h ls ) 

exp(-Phn) + exp(- ( 5/i i2 ) + exp(-/3/i i3 ) 

We minimize the following variational free energy (divided by the total number of spins N) 
with respect to h is : 

F = (<H > -TS )/N (6) 

The symbols < • • • > and S denotes respectively the expectation value and the entropy 
associated with the probability distribution P ({si}): 

<H> =J ]T E^( s )Pi( s )-7^ E Y,Pi( s )Pj( s ) ( 7 ) 

<i,j> ■n n s <i,j> nnn a 

^o = -EE^( s ) ln ^( s ) (8) 

i s 

Minimization problem with respect to h is is equivalent to that with respect to Pi(s) under 
following constraints: 

p i (l)+p j (2)+p i (3) = l , Pi {s)>0 (9) 
Translational symmetry assures that the free energy (0) is minimized when: 
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Pi ( s ) = C sA for all i 6 A (10) 
Pi ( s ) = C sB for all z E B (11) 

Obviously, C sa coincides with the expectation value of the concentration of s'th spin state 
on the sublattice a, calculated with the probability distribution Po{si}. They are under 
following constraints: 

Cla + C*2a + C*3a = 1 , C sa > (12) 

Then < H > and So are expressed with C sa : 

< H > = ^ £ C sA C sB - + C s 2 B )/2 (13) 

So = -N^Csa In C sA + C sB In C sB )/2 (14) 

where zx,z 2 denotes the coordination number of n.n. and n.n.n. sites respectively (for the 
simple cubic lattice, z\ = 6 and z 2 = 12). 

Furthermore, we define two-component sublattice magnetization similar as Ono used 



x a = (Ci a - C 2a )/V3 (15) 
y a = (C la + C 2a -2C 3a )/3 (16) 

Then the three quantities Ci a ,C2 a , and C 3a can be expressed by two quantities x a and y a 
owing to the constraint C\ a + C 2a + C 3a = 1. 

1 1 

Cia = - + —x a + -y a (17) 
1 v^3 1 

° 2a = 3 ~ !T Xq + 2^ a ( ' 18 ' ) 

C 3Q = i-^ (19) 

The two-component sublattice magnetization (x a , y a ) carry an irreducible, unitary represen- 
tation of the permutation group of the three spin states {1, 2, 3}. Owing to the constraints 



C sa > 0, the sublattice magnetization (x a ,y a ) is restricted to take a value within a regular 
triangle on the x a — y a plane (Fig. 2). Three vertexes of the triangle correspond to completely 
ordered state {C\ a = 1,C*2 Q = 0,Cs a = etc.) and the point (0, 0) corresponds to completely 
disordered state (C\ a = C*2 a = C% a = 1/3). 

We have minimized the free energy F with respect to xa,Va,%b, and %/b numerically 
using the following iteration: 



x (n + l) =x (n)_ A ^0 



a a 



(") „ _v< n ) 



dx Xa — ^ a — ^ c 



y(n+l) _ y(n) _ A ° | (on) 



where A is a small, positive quantity. The iteration fl2"0| ) is repeated until F converge to 
some minima. Several values were used as an initial value X^°) and to find absolute 
minimum value of Fq out of all minima. 

Finally we obtained a phase diagram of two parameters T j J and 7 (Fig. (3)). The order- 
disorder transition is continuous for 7 < 3/2 and discontinuous for 7 > 3/2. There are two 
kind of ordered phases, BSS phase and PSS phase, separated by a discontinuous transition 
line. Typical values of the concentrations of the three states and the sublattice magnetiza- 
tions of each ordered phases are shown in Fig. (4). In the PSS phase, the angle between each 
sublattice magnetizations is non-trivial value; it is greater than 120° and less than 180°. The 
non-trivial value of the angle can be understood as follows. The antiferromagnetic n.n. in- 
teraction makes two points (xa,Ua) and (xb,Vb) "repulsive", therefore prefers 180°. On the 
other hand, the energy of the ferromagnetic n.n.n. interaction is minimized when (x a ,y a ) 
locates on the vertex of the triangle, therefore prefers 120°. Thus we can understand the 
non-trivial value of the angle as a result of the competition of the two interactions. 

Banavar and Wu have studied the same model as @ with q = 3,4 using the mean-field 
theory [^TJ and concluded that the PSS phase does not appear. Our result disagree with 



theirs. 
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III. EFFECTIVE FREE ENERGY FORM 



Since the numerical method ( P0|) becomes less precise near the critical line, we expand 
the mean-field free energy F in powers of the order-parameter to investigate the critical 
behavior. 

Firstly we define ferromagnetic and antiferromagnetic order-parameters: 

X F = Xa + Xb , %AF = Xa~ %b 

Vf = Ua ~\~ Ub i Uaf = Ua Ub (21) 

Relevant quantities to the order- disorder transition are xaf and tjaf, which carry an irre- 
ducible representation of a group P3 (permutation of three spin states) xP 2 (permutation 
of two sublattices), isomorphic to the group c 6 „ (point group of a regular hexagon). Indeed, 
the allowed range of the antiferromagnetic order-parameter (xaf, Uaf) is a regular hexagon 
in the xaf — Uaf plane. 

Then we introduce polar coordinates of ferromagnetic and antiferromagnetic order- 
parameter: 

xp = R F cos9 F , y F = R F sm9 F 
xaf = Raf cos 9 A f , Vf = Raf sin 9 A f (22) 

The PSS phase and the BSS phase can be distinguished by the direction of the antiferro- 
magnetic order-parameter 9af', the value of 9af expected in the PSS and the BSS phase 
is kn/3 and {k + l/2)7r/3 {k = 0,1,2,3,4,5), respectively. Thus a quantity cos6#af is a 
relevant quantity to the PSS-BSS phase transition |14j]. The value of cos 69 af is —1 in the 
BSS phase, while cos 69 'af = 1 in the PSS phase. 

Now we use Raf,9af,Rf, and 9 F as independent variables of the free energy F and 
trace out R F and 9 F in order to obtain an effective free energy form which is expressed by 
Raf and 9af only. 

Faf(Raf, 9af) = mm F (Raf, 9af, Rf, 9 f ) = F (Raf, 9af, Rf, 9 f ) (23) 



where Rf and Of gives minimum value of Fq for fixed Raf and Oaf- Location of minima 
Rf and Of is determined by solving the following equations: 

dF , 



q r \{R F ,e F )=(R F ,e F ) 



dF 0l 



(24) 



(25) 



QQ F \{R F ,e F )=(R F ,e F ) 

Since we cannot solve ([24]) and (|25|) explicitly, we expand -Rf and $f in powers of Raf which 
is small around the critical line. Firstly we expand ( |24] ) and (^) in powers of both Rp and 
Raf- Lowest order terms follow: 

dF _ 9(2 J - 4 7 J + T) 27T sm(20 AF + B F ) d2 

dR~ F ~ 8 64 AF 

+0{R 2 F ) + 0(i^ F ) + 0{RfR 2 af ) (26) 

<9F / 27Tcos(20 AF + F ) 2 \ 
9$f ^ 64 Raf ) 

+R F (0(R 4 af ) + 0(R%) + 0{R\ F R F j) (27) 

Equations (|26|) and fl27l) indicate that i? F ~ 0(i?^ F ) and cos(26*af + #f) ~ C(i?^ F ), so 
we assume that Rp and Op can be expanded as below: 

Rp = c%R 2 af + c aRaf + c 6-R^f ' ' ' (28) 
F = vr/2 - 20 AF + d 2 R 2 AF + d 4 R AF + d 6 R AF ■ ■ ■ (29) 

The coefficients can be determined by letting (|28|) and (|29|) into equations (^4j) and (^5|), 
leading to the following result: 

c 2 = 3r/8 (30) 

d 2 = A(_ 2 + 4r-r 2 )sin(6^F) (31) 

C4 = g(r - r 2 ) + ^(-2r + 4r 2 - r 3 ) cos(60 A f) (32) 
rf 4= 2 ^(-l + 4r-4r 2 + r 3 )sin(60 A F) 

+ g i §(2-4r + r 2 ) 2 sin(120 A F) (33) 
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where t = T/J and r = t/(2 - 47 + t). 



Finally the effective free energy form is obtained by letting (|28| ) and ( f29"|) into (|23|) . Note 
that only the terms which is invariant under the transformations of Cq v are present in Faf, 
and they can be written as R A n £ 2m cos{6n9 AF ) , (n, m > 0). We have calculated the effective 
free energy F AF up to sixth order term as below: 

F AF = A + A 2 R\ f + A 4 R AF + (A 6 + B 6 cos 69 AF )R 6 AF + 0{R 8 AF ) (34) 

where 

A = J-2 7 J-Tlog3 (35) 

A 2 = JLj(_ 2 -4 7 + t) (36) 

A, = T | i T(2 - r) (37) 

A 6 = f|T(4 - 6r + 3r 2 ) (38) 

fi 6 = ifloTM + 30r - 30r2 + 5r3 ) ( 39 ) 

The critical line is the region A 2 = and A4 > 0, which corresponds to the region 
t = 2 + 47, 7 < 3/2. The order-disorder transition is discontinuous for 7 > 3/2 and 
(7,t) = (3/2,8) is a tricritical point where A 2 and A4 vanish simultaneously |22] . 



The six-fold anisotropy term BqR af cos is the origin of the BSS and the PSS phase; 
positive B 6 corresponds to the BSS phase and negative B 6 to the PSS phase. Presence of 
higher order anisotropic terms such as R\ 2 f cos ^-2@af allows a occurrence of minima of Faf 
at Oaf 7^ kn/6, which corresponds to neither the PSS nor the BSS phase. However, the 
positions of minima remain at 9af — kn/6 if the coefficients of such higher order anisotropic 
terms are sufficiently small. The following example helps understanding above issue: a 
positions of minima and maxima of a function f(6) = cos 6 + a cos 26 remain at 9 = kn, as 
long as \a\ < 1/4. 

On the critical line t = 2 + 47, B 6 is expressed as a function of 7: 

B 6 = 243(2 + 4 7 ) (1 + 3O7 - I8O7 2 + 407 3 )/1310720 (40) 



It vanishes at 7 ~ 0.2032 on the critical line, where the BSS-PSS transition line merges into 
the critical line. 

Note that higher order term is needed to investigate off-critical regions, where the nu- 
merical method ( |20|) works well. So we do not calculate higher order terms. 



IV. CONCLUSION 

In this paper we have shown that a mean-field theory analysis indicates the presence of 
a tricritical point in the three-state Potts model with the antiferromagnetic n.n. interaction 
and the ferromagnetic n.n.n. one. It is believed that mean-field theory is qualitatively 
correct above the upper critical dimension, where spatial fluctuations of order-parameter 
at the critical region become negligible. Since the upper critical dimension of the tricritical 
phenomena is three ||22|| , we expect the model on the simple cubic lattice to posses a tricritical 
point. Owing to the simpleness of the Hamiltonian, this model must serve to develop 
numerical methods for studies of a tricritical phenomena, such as finite-size scaling in Monte 
Carlo simulations. 

As for the nature of the ordered phase, we have shown that the PSS phase is realized 
at the strong n.n.n. coupling region, as a result of the competition of the two kind of 
interactions, while the BSS phase is realized at the weak n.n.n. coupling region. 

However, it should be noted that mean-field type analysis like this work may not give 



collect information, as pointed in Ref. ||14j| , about the six- fold anisotropy term which becomes 
relatively small compared to the order-parameter fluctuation near the critical line. So further 
study, such as Monte Carlo simulations, may be needed to clarify the nature of the ordered 
phase just below the critical line. 
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Figure Captions 

Fig.l. Typical ground state configuration of antiferromagnetic three-state Potts model on 
a square lattice. Two different symbols ( Dand o) show the two interpenetrating sublattices. 

Fig. 2. A range of sublattice magnetization (x a , y a ). Three vertexes 1,2, and 3 correspond 
to the completely ordered states in which all the spins in the sublattice take the same value 
1,2, and 3, respectively. 

Fig. 3. Schematic phase diagram of two parameters T/J and 7. The solid line indicates 
a continuous transition and the dashed line indicates a discontinuous one. The continuous 
order-disorder transition line ends at 7 = 1.5. The BSS-PSS phase boundary merges into 
the order-disorder transition line at 7 ~ 0.2032. 

Fig. 4. Typical values of the concentrations of the three states (upper) and the sublattice 
magnetizations (lower) of each ordered phases. 
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Fig. 2. A range of sublattice magnetization (x a} y a ). Three vertexes 1,2, 
and 3 correspond to the completely ordered states in which all the spins in 
the sublattice take the same value 1,2, and 3, respectively. 
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Fig. 3. Schematic phase diagram of two parameters T j J and 7. The 
solid line indicates a continuous transition and the dashed line indicates a 
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1.5. The BSS-PSS phase boundary merges into the order-disorder transition 
line at 7 ~ 0.2032. 
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Fig. 4. Typical values of the concentrations of the three states (upper) 
and the sublattice magnetizations (lower) of each ordered phases. 
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